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Abstract 

Consider the wave equation with constant or variable coefficients in IR 3 . The ini- 
tial datum is a random function with a finite mean density of energy that also satisfies 
a Rosenblatt- or Ibragimov-Linnik-type mixing condition. The random function con- 
verges to different space- homogeneous processes as X3 — > ±00 , with the distributions 
fj,± . We study the distribution at of the random solution at a time t G IR . The main 
result is the convergence of fit to a Gaussian translation-invariant measure as t — > 00 
that means central limit theorem for the wave equation. The proof is based on the 
Bernstein 'room-corridor' argument. The application to the case of the Gibbs mea- 
sures fj,± = g± with two different temperatures T± is given. Limiting mean energy 
current density formally is —00 • (0, 0, T + — T_) for the Gibbs measures, and it is finite 
and equals to — C(0, 0, T + — XL) with C > for the convolution with a nontrivial test 
function. 
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1 Introduction 



The paper concerns a mathematical problem of foundations of statistical physics. We con- 
sider the Second Law of thermodynamics in a reversible infinite dimensional Hamiltonian 
equations. The Law states that the energy current is directed from higher temperature to 
lower temperature. We derive the Law for wave equations in IR 3 with constant and variable 
coefficients. The key role plays the mixing condition of Rosenblatt- or Ibragimov-Linnik- 
type for an initial measure. The mixing condition is introduced initially by R.L. Dobrushin 
and Ya.M. Suhov in their approach to the problem of foundation of statistical physics for 
infinite-particle systems, [4, 5]. The mixing condition is used also in the paper [2] which con- 
cerns a discrete version of our result for a ID chain of harmonic oscillators. Let us explain 
our result in the case of constant coefficients, 



u(x, t) = Au(x, t), x e IR , 



u 



t=0 



u (x), u 



(1.1) 



= v {x) 

t=0 y ' 



Denote Y(t) = (Y°(t), Y 1 ^)) = (u(-,t),u(;t)) , Y = (Y °X) = (u ,v ) . Then (1.1) 
becomes 

Y(t)=F(Y(t)), teJR, Y(0) = Y . (1.2) 

We assume that the initial datum Yq is a random function with zero mean living in a 
functional phase space H of states of finite local energy; the distribution of Y is denoted 
by fio . Denote by ji t (dY) , t G IR, the measure on H giving the distribution of the 
random solution Y(t) to problem (1.2). We assume that the initial correlation functions 
Q l J(x,y) = e(Yq{x)Yq (y)^j , i,j = 0,1, and some of their derivatives are continuous and 
decaying as \x — y\ — > oo . In particular, the initial mean energy density is bounded: 

£[|VM (a:)r + K(x)| 2 ] = [V* • V y Q° \x, y)\\ v=x + QR{x,x) < C < oo, x G IR 3 . (1.3) 
Next, we assume that the initial correlation matrix (Qq (x, y))jj=o,i has the form 

f ql(x-y), x 3 ,y 3 <-a, 
Ql 1 (x,y) = l (1.4) 
[ q+(x-y), x 3 ,y 3 > a. 

Here q±(x — y) are the correlation functions of some translation-invariant measures fi± with 
zero mean value in Ti , x — (x±, X2, x 3 ), y = (1/1,1/2,1/3) G IR 3 , and a > 0. The measure 
Ho is not translation-invariant if q l l 7^ q % l . Finally, we assume that the initial measure /xo 
satisfies a mixing condition. Roughly speaking, it means that 

Yq(x) and Y (y) are asymptotically independent as \x — y\ ^ 00. (1.5) 

Our main result establishes the (weak) convergence 

Ht^Hoo, t^oo, (1.6) 

to an equilibrium measure /ioo , which is a translation-invariant Gaussian measure on H . 
The similar convergence holds for t — > — 00 since our system is time-reversible. We construct 
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generic examples of the random initial datum satisfying all assumptions imposed. We get 
the explicit formulas (2.13)-(2.15) for the limiting correlation matrices. 

We apply our results to the case of the Gibbs measures /i± — g±. Formally 

1 /(|Vtx (a;)| 2 + |uo(a:)| 2 )<fo TT 

g±(du ,dv ) = —e 2 J \ [du (x)dv (x), f3± = T± , (1.7) 

Z ± X 

where T± > are the corresponding absolute temperatures. We adjust the definition of 
the Gibbs measures g± in Section 3. The Gibbs measures g± have singular correlation 
functions and do not satisfy our assumptions (2.10). Respectively, our results can not be 
applied directly to g± . We reduce the problem by a convolution with a smooth function 
9 E D = C^°(IR 3 ) : we consider Gaussian processes u± corresponding to the measures g± 
and define the "smoothened" measures g± as the distributions of the convolutions u± * 9 . 
The measures g± satisfy all our assumptions, and the convergence gf —r g 6 ^ follows from 
(1.6). This implies the weak convergence of the measures g t —r g^ since 9 is arbitrary. We 
show that the limit energy current for g^ is formally 

J oo = -oo-(0,0,T + -T_). 

The infinity means the "ultraviolet divergence". This relation is meaningful in the case of 
smoothened measures g^ , 

JL = -CV(0,0,T + -T_), 

if 9{x) is axially symmetric with respect to Ox 3 ; Cq > if 9{x) ^ . This corresponds to 
the Second Law of thermodynamics. 

We prove the convergence (1.6) in three steps using the strategy of [10, 18, 19]. 

I. The family of measures /i t , t > , is weakly compact in an appropriate Frechet space. 

II. The correlation functions converge to a limit: for i,j = 0, 1, 

y) = J Y i {x)Y\y) l x t {dY) - QZ(x,y), t - oo. (1.8) 

III. The characteristic functionals converge to the Gaussian: 

= JeMi(Y,*)ht(dY)^eM-lQoo(y,*)}, t ^ oo, (1.9) 

where \& is an arbitrary element of the dual space, and Qoo is the quadratic form with the 
integral kernel (QZ(x, y))ij=o,i ■ 

Property I follows from the Prokhorov Compactness Theorem by using methods of [23]. 
First, one proves a uniform bound for the mean local energy with respect to the measure 
fit ■ The conditions of Prokhorov' Theorem then follow from Sobolev's Embedding Theorem. 
We deduce the uniform bound from the explicit expression for the correlation functions 
Ql J (x,y) . The expression follows from the Kirchhoff formula for the solutions to (1.1). In 
particular, in the case u (x) = , we have 

u(x,t) = ^- t J v (x')dS(x'), (1.10) 

Si Or) 
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where dS(x') is the Lebesgue measure on the sphere S t (x) : \x' — x\ — t . 

Property II also follows from explicit formulas for Ql 3 (x,y) . The formula (1.10) allows 
to express the correlation functions Q l t J (x, y) in terms of integrals over spheres of radius t . 
In the limit, t — > oo , the spheres become the planes. Respectively, Q t ^ {x 1 y) is expressed 
in terms of integrals of the Radon transform of initial correlation functions Qq (x, y) . We 
reduce the expressions to some convolutions. 

Remarks i) The dynamics (1.1) is translation invariant, and its Fourier transform has a 
very simple form. However, the proof of (1.8) in Fourier transform is not transparent and 
requires additional efforts since our main assumption (1.4) is stated in the coordinate space, 
ii) Our proof of the convergence (1.8) in Sections 5 and 6 does not allow a simplification in 
the particular case of the Gibbs measures (1.7). This is related to the slow long-range decay 
of the correlation function Qq°(x, y) ~ \x — , \x — y\ ^ oo . 

We deduce property III using the method of [10]. The method is based on a modi- 
fication of the Bernstein room-corridor' argument, and it is suggested by the structure of 
the Kirchhoff formula (1.10): roughly speaking, (1.10) is "the sum" of weakly dependent 
random values devided by the square root of their "number" . This observation allows us to 
reduce the proof of (1.9) to the Lindeberg Central Limit Theorem, similarly to [10]. We do 
not consider the case n = 2 : it requires a different approach since the strong Huyghen's 
principle breaks down. 

Let us note that our mixing condition is weaker than that in [10]: this is necessary in the 
application to the Gibbs measures (1.7). Namely, we introduce different mixing coefficients 
for partial derivatives of the random solution at t — : we assume that the long-range decay 
of the mixing coefficients depends on the order of the derivatives. Respectively, our proof 
requires new tools (see Sections 7, 9, 10). For instance, the splitting (7.15) and the bound 
(9.5) play a crucial role. 

All the three steps I-III of the argument rely on the mixing condition. Simple examples 
show that the convergence to a Gaussian measure may fail when the the mixing condition 
fails (see [10]). 

In conclusion, we extend the convergence in (1.6) to the equations with variable coeffi- 
cients, that are constant outside a finite region. The extension follows immediately from our 
result for constant coefficients, using method of [10]. The method is based on the scattering 
theory for the solutions of infinite global energy, which is constructed in [10]. 

The paper is organized as follows. In Section 2 we formally state our main result. We 
apply it to the Gibbs measure in Section 3. Sections 4-10 deal with the case of constant coef- 
ficients: the compactness (Property I) and the convergence (1.8) are proved in Sections 4-6. 
In Section 7 we introduce the 'room-corridor' method, in Section 8 we prove the convergence 
(1.9), and in Sections 9, 10 we check the Lindeberg condition. In Section 11 we establish 
the convergence (1.6) for variable coefficients. Appendix A concerns the Radon transform 
and convolution, and Appendix B concerns the Gaussian measures in the weighted Sobolev 
spaces. 

Let us note that the equation (1.1) describes a continuous n -dimensional family of har- 
monic oscillators. Therefore, our result is an extension of the results [2, 21] that concern the 
infinite one-dimensional chains of harmonic oscillators. 
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Our formulas for the limit correlation functions correspond to the discrete one-dimensional 
version [2]. For instance, the position- momentum correlations have a power long-range de- 
cay. On the other hand, in [20, 16] the limit correlation functions are constructed for the 
finite chains of N oscillators with the "Langevin" boundary value conditions. In the limit 
iV — + oo the correlation functions have an exponential long-range decay. This means that 
this limit leads to another stationary measure of the infinite chain, different from [2, 21]. 

The convergence to statistical equilibrium for the wave equation is established in [10] (see 
also [18, 19]) for the case of a translation-invariant initial measure /x • This corresponds to 
our result in the particular case T_ = T + . The similar result has been proved for the Klein- 
Gordon equation, [9, 14]. If the initial measure /j,q coincides with one of the equilibrium 
limit measures //oo , the corresponding random solution Y(t) is mixing in time, [6, 7, 8]. 

2 Main results 
2.1 Notations 

We assume that the initial datum Y belongs to the phase space 7i defined below. 

Definition 2.1 H = #/ oc (IR 3 ) © if,° oc (IR 3 ) is the Frechet space of pairs Y = (u(x),v(x)) 
of real functions u(x) , v(x) , endowed with the local energy seminorms 

\\Y\\ 2 R = J {\u{x)\ 2 + \Vu{x)\ 2 + \v{x)\ 2 )dx< oo, Vi?>0. (2.1) 

\x\<R 

Proposition 2.2 follows from [15, Thms V.3.1, V.3.2] as the speed of propagation for Eqn 
(1.1) is finite. 

Proposition 2.2 i) For any Y G 7i there exists a unique solution Y(t) G C(IR, Ti.) to the 
Cauchy problem (1-2). 

ii) For any t G IR, the operator U(t) : Y i— > Y(t) is continuous in Ti . 
Hi) The energy inequalities hold \/R > , 

||^)y ||fl<C(*)||yo||fl+| t |, teJR. (2.2) 

Let us choose a function ((x) G C^°(1R 3 ) with C(0) ^ 0. Denote by Hf oc (M 3 ), s G IR, 
the local Sobolev spaces, i.e. the Frechet spaces of distributions u G //(IR 3 ) with the finite 
seminorms 

\\u\\ s ,r := \\A s (({x/R)uj\\ L 2 {R3) , 

where A s v := F^ x ((k) s v(k)) , (k) := yj\kfi + 1 , and v := Fv is the Fourier transform of a 
tempered distribution v . For G D = C^°(IR 3 ) define Fip(k) = J e ik ' x ^{x)dx. 

Definition 2.3 For s G IR denote W = ^"(IR 3 ) © #f oc (IR 3 ). 

Using the standard techniques of pseudodifferential operators and Sobolev's Theorem 
(see, e.g. [13]), it is possible to prove that 7i° = 7i C 7i~ £ for every e > 0, and the 
embedding is compact. We denote by (•, •) the scalar product in real Hilbert space L 2 (IR 3 ) 
or in L 2 (IR 3 ) © IR^ or in its various extensions. 
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2.2 Random solution. Convergence to equilibrium 

Let (Q, a, P) be a probability space with the expectation E , let B(H) denotes the Borel a - 
algebra in H . We assume that Y = Y (ui, x) in (1.2) is a measurable random function with 
values in (H, B(H)) . In other words, (cu,x) i— > Y (cu,x) is a measurable map n x n 3 -> ir 2 
with respect to the (completed) ct -algebra E x £(IR 3 ) and £(IR 2 ) . Then F(t) = E/"(*)y 
is also a measurable random function with values in (7i,£>(7i)), due to Proposition 2.2. 
We denote by fio(dY ) the Borel probability measure in 7i that is the distribution of Y . 
Without loss of generality we assume (Q, E, P) = {7i, B(H), /io) and Y (u,x) = u(x) for 
Ho(du) x dx -almost all (oo,x) G H x IR 3 . 

Definition 2.4 /x t zs £/ie Borel probability measure in 7i that is the distribution of Y(t) : 

^(B) = fio(U(-t)B), VB G BiU), t G H. (2.3) 

Our main goal is to derive the convergence of the measures /x t as t — > oo . We establish 
the weak convergence of /j, t in the Frechet spaces Hr e with any e > : 

A*t — ^ ^ as t — > oo, (2.4) 

where /Xoo is the Borel probability measure in the space . By definition, this means the 
convergence 

J f(Y)nt(dY) - | /(y) Aloo (dy) as t - oo (2.5) 
for any bounded continuous functional f(Y) in the space H~ £ . 

Definition 2.5 The correlation functions of the measure /x t are defined by 

Q i t j (x,y) = E(Y i (x,t)Y j (y,t)), i,j = 0,1, for almost all x,?/GlR 3 xIR 3 (2.6) 
if the expectations in the RHS are finite. 

We set V = D®D, and (Y, = (Y°, *°) + (y 1 ,^ 1 ) for y = (y ^ 1 ) G H, and 
\j/ — (\j/0 ) xp 1 ) g x> . For a Borel probability measure /x in the space 7i we denote by fx the 
characteristic functional (Fourier transform) 

/}(*) = /" exp(i(y *)) /i(rfy), $GD. 

A measure /x is called Gaussian (with zero expectation) if its characteristic functional has 
the form 

/}(*) = exp{-^Q(*,*)}, *GD, 

where Q is a real nonnegative quadratic form in V . A measure /x is called translation- 
invariant if 

n(T h B) = n(B), V5 G /i G H 3 , 

where T ft y(x) = Y(x — h) . 
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2.3 Mixing condition 

Let 0(r) denote the set of all pairs of open subsets A, B C 1R 3 of distance p(A, B) > r , 
let a = («!, ct2, a 3 ) with integers ctj > . Denote by a ia (A) the a -algebra of the subsets 
in TL generated by all linear functionals 

Yt-+(D a Y i ,ip) = J D a Y i (x)ip(x)dx, |a|<l-i, i = 0,l, 



JR: 
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where ip G D with supp?/> C .4.. For d — 0, 1 let be the a -algebra generated by o" ia 
with i + |a| > d , i.e. 

o"d = V a ia, d = 0,l. 

i+\a\>d 

We define the Ibragimov-Linnik mixing coefficient of a probability measure /io on H. (cf. 
[11, Dfn 17.2.2]) for d u d 2 = 0, 1 as 

0d lld2 (r)= sup sup — . 

> 



Definition 2.6 TTie measure /io satisfies the strong uniform Ibragimov-Linnik mixing con- 
dition if for any di,d 2 = 0, 1 

<f> dl , d2 (r)^0, r^oo. (2.7) 
Below we specify the rate of the decay. 

2.4 Main theorem 

Let Vd £ C[0, oo) denote some continuous nonnegative nonincreasing functions in [0, oo) 
( d — 0, 1, 2 ) with the finite integrals, 

oo 

J (1 + r) d "V d (r)dr < oo. (2.8) 
o 

We also denote v(r) = v 2 (r) . We assume that the measure /i satisfies the following 
conditions S0-S3: 

50 /io has the zero expectation value, 

EFoOr) = 0, x G 1R 3 . (2.9) 

51 The correlation functions of /i have the form (1.4). 

52 The following derivatives are continuous and the bounds hold, 

mm*, »>i<{ £<j* : °<± ; | <-. + , +w+ ia. (2.io, 

53 The measure /io satisfies the strong uniform Ibragimov-Linnik mixing condition, and for 
d±,d 2 = 0, 1 

dl ,d 2 (r) <CV d 2 (r), d = d! + d 2 . (2.11) 
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Remark 2.7 i) Condition S2 implies (1.3). Condition S3 implies the estimates (2.10) with 
i + \a\<l, j + \P\<l. 

ii) The conditions S2 and S3 allow various modifications. We choose the variant which 
allow an application to the case of the Gibbs measures (1.7) (see the next section). Our 
mixing condition S3 is weaker than the mixing condition [10] which corresponds to S3 with 
the functions ^o,i(r) < v i{ r ) ■ On the other hand, the estimates (2.10) with d > 2 are not 
required in [10]. 

Let Six) = j — | be the fundamental solution of the Laplacian, i.e. AS = Six) for 

Att\x\ 

x E IR 3 , and P(x) = —iF^ 1 — — — where F' 1 is the inverse Fourier transform. Define, for 

| A; | 

almost all x, y G IR 3 , the matrix- valued function 

Qoo(x,y) = (Q* > (x,y)).. = - y)) . . , (2.12) 

V / i,j=0,l \ /jj=0,1 



where 



oo 



\ [ql° + q 00 -S* + q]}) + P * - q 01 - q™ + q W )} , (2.13) 
1 r ,10 , U0 „01 „01 



q ™ + gi" _ q f _ ^ + p * _ g ii _ Agf + Ag uu )j , (2.14) 

<£ = -A<£ = ^ [tf + d 1 - AGff + g 00 ) + P*A(gf-gi°-g° 1 + g 01 )]. (2.15) 

The definition of the convolutions with P in formulas (2.13)-(2.15) is adjusted in Appendix 
A (formula (6.10)). 

Denote by Qoo(^, the real quadratic form in V defined by 

Q O0 {*,*) = E / Ql(x,yW(x)V(y)dxdy. (2.16) 

ij=0,1 R 3 xR 3 

Our main result is the following theorem. 

Theorem 2.8 Let S0-S3 hold. Then 

i) the convergence in (2.4) holds for any e > . 

ii) The limiting measure fi^ is a Gaussian equilibrium measure on 7i . 
Hi) The limiting characteristic functional has the form 

= ex P {-is oo (^, *)}, 

where <2oo is the quadratic form with the integral kernel Qoo(x, y) defined in (2.12)-(2.15). 

Theorem 2.8 can be deduced from Propositions 2.9 and 2.10 below, by using the same 
arguments as in [23, Thm XII. 5. 2]. 

Proposition 2.9 The family of the measures {/i t , t > 0} , is weakly compact in Hr £ with 
any e > . 

Proposition 2.10 For any $6P, 

£*(*) = J exp(t(Y,^)) fi t (dY) ^ expi-^Q^,^)}, t - oo. (2.17) 

Proposition 2.9 is proved in Section 4 for a simple particular case, and in Section 6 for 
the general case. Proposition 2.10 is proved Sections 7, 8. 
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2.5 Examples 



2.5.1 Gaussian measures 



We construct the Gaussian initial measures /xo satisfying SO S3. Let us take some Gaussian 
measures /x± in Ti with correlation functions q± (x — y) which are zero for i ^ j , while for 

2 = 0,1, 



'l + \k\) s d2ql(k) e L^JR 3 ), < d = 2i + s < 4, | 7 |<l + d, 
qf(k) > 0. 



(2.18) 



Then fj,± satisfy SO, S2 with the functions va(r) — C(l + r) 1 d for a sufficiently large 
C > . Let us take the functions C± e C°°(IR) s.t. 



Us) 



1, for ± s > a, 
0, for ± s < —a. 



Let us introduce (YL, i+) as a unit random function in the probability space (7i x Ti, /j_ x 
//+) . Then Y± are Gaussian independent vectors in . Define /io as the distribution of 
the random function 

Y (x) = C-MY-(x) + C+(x 3 )Y + (x). (2.19) 
Then correlation functions of /jq are 



Qo( x ,y) = Q-( x ~ y)C-MC-(y 3 ) + q+(x - y)(+(x 3 )( + (y 3 ), i,j = 0,1, 



(2.20) 



where x = (xi,X2,x 3 ) , y = (1/1,1/2,1/3) £ 1R 3 , 9± are the correlation functions of the 
measures /i± . Then SO and SI hold, and S2 follows for /x with the same functions z/rf(r) 
as for /j,± . Let us assume, in addition to (2.18), that 



q±(x) = 0, \x\ > r . 



[2.21) 



Then the mixing condition (2.7) holds since 4>d 1 ,d 2 ( r ) = 0, r > r , and S3 follows. For 
instance, (2.18) and (2.21) hold if qt(h,k 2 ,k 3 ) = f(k 1 )f{k 2 )f{k 3 ) with 

f(z) = ((1 - cos(r z/VS))/z 2 ) N , zeK, 
where iV > is an integer, 2iV — s > 1 ( s = 4 — 2i ). 



2.5.2 Non-Gaussian measures 

Let us choose some odd nonconstant functions /°, J 1 G C 4 (IR) with bounded derivatives. 
Let us define /Xq as the distribution of the random function (/°(F°(a;)), / 1 (F 1 (a;))), where 
(yo^y 1 ) is a random function with the Gaussian distribution /x from the previous example. 
Then SO, SI and S3 hold for /Xq with some appropriate functions ^ since corresponding 
mixing coefficients <p* dl d2 (r) = for r > r . Therefore, SO implies for the corresponding 
correlation functions Q^(x,y) = for \x — y\ > r , so S2 also holds. The measure /Xq is 
not Gaussian since the functions f° , f 1 are bounded and nonconstant. 



8 



3 Application to Gibbs measures 



We apply Theorem 2.8 to the case when /i± are the Gibbs measures (1.7) corresponding to 
different positive temperatures T_ ^ T + . 

3.1 Gibbs measures 

We will define the Gibbs measures g± as the Gaussian measures with the correlation func- 
tions (cf. (1.7)) 

q™(x-y) = -T ± £(x-y), q£(x - y) = T ± S(x - y), q<£(x - y) = q™(x - y) = 0, (3.1) 

where x, y G H 3 . The correlation functions q± do not satisfy condition S2 because of 
singularity at x = y . The singularity means that the measures g± are not concentrated in 
the space Ti . Let us introduce appropriate functional spaces for measures g± . First, let us 
define the weighted Sobolev space with any s,a6E. 

Definition 3.1 H Sta (JR 3 ) is the Hilbert space of the distributions u G S"(IR 3 ) with the finite 
norm 

\\u\\ 8 , a = || (x) Q A a «|| £a(Il 3 ) < oo, A s u = F- 1 [(k) s u(k)] . (3.2) 
Let us fix arbitrary s, a < — 3/2 . 
Definition 3.2 G s>a is the Hilbert space iJ s+liQ ,(IR 3 ) © H s ^ a (1R 3 ) , with the norm 

Wts,a = \\ u \\s+l,a + \\v\\ 8 ,a < OO, Y = (u,v). 

Introduce the Gaussian Borel probability measures g±(du) , g±(dv) in spaces H s+ i :a (]R 3 ) 
and H s>a (]R 3 ) , respectively, with characteristic functionals 



0±(VO = J exp{z(u, ij>)}g±(du) = exp{ ^ 2 ^ ^ } 
9±W = J exp{i(v,i>)}gl(dv) = exp{— ^|^} 



i> G D. 



By the Minlos theorem, [3], the Borel probability measures g± , g± exist in the spaces 
H s+ i :(x (M 3 ) , H S:a (M 3 ) , respectively, because formally (see Appendix B) 

/ \\ u \\ 2 s+i, a gl( du ) < oo, J \\v\\l a gl(dv) < oo, s,a < -3/2. (3.3) 

Finally, we define the Gibbs measures g±{dY) as the Borel probability measures g°±(du) x 
g±{dv) in G s ^ a . Let g^idY) be the Borel probability measure in G s>a that is constructed 
as in the Example of previous section with /j,±(dY) = g±(dY) . It satisfies SO and SI with 
q± from (3.1). However, g does not satisfy S2. Therefore, Theorem 2.8 cannot be applied 
directly to Ho = go ■ Gs,a C H s by the standard arguments of pseudodifferential equations, 
[13]. The next lemma follows by Fourier transform from the finite speed of propagation for 
wave equation. 

Lemma 3.3 The operators U{t) : Y \— > Y(t) allow a continuous extension H s i— > H s . 
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3.2 Convergence to equilibrium 



Let Y be the random function with the distribution g , hence Y e G s ^ a a.s. Denote by 
g t the distribution of U(t)Y . 

Theorem 3.4 Let s < — 5/2 . Then there exists a Gaussian Borel probability measure 
in 7i s such that 

9t — 7 goo, t^oo. (3.4) 

Proof Let us fix an s < — 5/2 and introduce the random function F S := A S Y , F S e Go,a 
a.s. Let us denote by gf the distribution of U(t)(A s Y ), t G IR. Then <jf = g t A~ s , 
and ^ = g*A s since A s (£/(t)F ) = U(t)(A s Y Q ). Let us denote by Q s t (x,y) the (matrix) 
correlation function of measure gl . 

Measure g$ obviously satisfies SO. The correlation function Q s Q (x,y) also satisfies SI 
with a suitable modification: (1.4) holds up to 5(1 + \x\ + \y\)~ N with any 5, N > and 
with a = a(<5) . This follows from the convolution representation Qq(x, y) = Qo(x,y) * 
(A s (x)A s (y)) since A s (x) = F" 1 (k) s is a function e Lj^IR 3 ) with a rapid long-range 
decay. S2 also holds for #q with the functions ^(r) = C(l + r) _1_(i for a sufficiently 
large C = C(s,T±) > . It follows immediately for s < with sufficiently large \s\ from 
the same convolution representation. For s < —5/2 it follows by the pseudodifferential 
operators techniques. 

Then the conclusions of Lemmas 4.1, 5.1 hold for the random function Yq and the 
correlation functions Qf(x,y) of the measures g% . The proofs are almost unchanged. Hence, 
the convergence (2.4) holds for the Gaussian measures gl : Ve > 

9 S t <4, * -> °o, (3-5) 
where g^ is a Gaussian measure in Ti . Therefore, 

9t — r goo, t-> oo, 

since gt = glA s . This implies Theorem 3.4. □ 

The limiting measure g^ is Gaussian with the correlation matrix = (Q^x, y))j J= o,i, 
where 

Q£(x,y)= <£(x-y)=-±(T + + T-)E(x-y), (3.6) 
Q™(x,y) = -Qfc(x,y)= tf»(x - y) = \{T + - T_)P{x - y), (3.7) 
Q^y)= q£(x-y)= ±(T + + T_)5(x-y). (3.8) 

The identities (3.6)-(3.8) follow formally from (3.1) and from (2.13)-(2.15). For the proof 
we apply (2.13)-(2.15) to the initial measure g$ . 
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3.3 Limit energy current density 

Let u(x,t) be the random solution to (1.1) with the initial measure /xq satisfying SO S3. 
The mean energy current density is Ej(x,t) = —Eu{x,t)\/u{x,t) . Therefore, in the limit 
t — > oo , 

£j(MWoo = v^(o). 

Respectively, in the case of the "Gibbs" initial measure go , the expression (3.7) for the 
limiting correlation function implies formally that 



T — T 

loo = VP(0), 



k sgn k 3 



where [VP] (2) = — F 1 — — — (z) . Hence, formally we have the "ultraviolet diverging" 
limit mean energy current density, 

T+ - T_ r k sgn k 3 



/^pdfc = -oo.(o,o,r + -r_). 



2(2tt) 3 

R 3 

On the other hand, for the convolution U (t) (Yq * 9) the corresponding limiting mean energy 
current density is finite, 

t = / IWI 2 ^* = -Co • (o,o,T + -r_), 

R 3 

if 9(x) is axially symmetric with respect to Ox 3 ; Co > if 9{x) ^ . 

4 Compactness of the measures family 

Proposition 2.9 can be deduced from the bound (4.1) below with the help of the Prokhorov 
Theorem [23, Lemma II. 3.1] as in [23, Theorem XII. 5. 2]. 

Lemma 4.1 Let S0-S2 hold. Then the following bounds hold 

svvE\\U{t)Y Q \\ 2 R < oo, R > 0. (4.1) 

Proof Assumption S2 and Proposition 2.2 iii) imply by the Fubini Theorem the exis- 
tence of the correlation functions in (2.6), where Y l (x,t) are the components of Y(x,t) = 
(Y°(x, t), Y 1 (x, t)) . Therefore, Definition (2.1) implies 



E\\Y{-,t)\\l = E J \Y°{x,t)\ 2 dx + E J \VY°{x,t)\ 2 dx + E J \Y\x,t)\ 2 dx 

\x\<R \x\<R \x\<R 

= J Q° t °(x,x)dx+ J V x -V y Q° t °(x,y)\ y=x dx+ J Q\\x,x)dx. (4.2) 

\x\<R \x\<R \x\<R 

We bound for example the integral of Q^°(x, x) in (4.2) in the particular case when Y Q ° = 
u (x) = almost surely. The general case will be considered in Section 6 as well as the 
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bounds for two remaining integrals in (4.2). Let us assume for a moment that the function 
Yq = vq is continuous almost surely. Then the Kirchhoff formula (1.10) gives by the Fubini 
Theorem, 

(ff( XyX ) = -±- J Ql\x'y)dS(x')dS(x"). (4.3) 

St(x)xS t (x) 

Let us assume for a moment that 

Qi j (x',x") = for \x'-x"\>r , i,j = 0, 1. (4.4) 
Then (4.3) implies the uniform bound 

Qt°(x,x) < ^ J dS(x')dS(x") <I = drf, t G JR. (4.5) 

S t (x) x S t (x) 
\x' — x"\ < r 

Hence, the bound follows, 

/ Q°°(x,x)dx <CIR 3 , teJR. (4.6) 

J\x\<R 

Next we remove the additional assumption (4.4) by the following known lemma on spherical 
integral identity, [12]. 

Lemma 4.2 Let h(r) G C(0, +oo). Then for any r > and x" G S t (x) the identity 
holds, 

2t 

J h(\x' - x"\)dS(x') = 2n J rh(r)dr. (4.7) 

{x'<=St(x): \x'-x"\>r } r 

Therefore, (4.3), S2 with d = 2 and Lemma 4.2 with r = imply (see (2.8)), 

Q™(x,x)<—— J v 2 (\x'-x"\)dS(x')dS(x")<C 1 J rv 2 (r)dr<C 2 <oo. (4.8) 

71 St(x)xS t (x) 

Then (4.6) follows without the assumption (4.4). The assumption on the a.s. continuity of 
v (x) can be removed by a convolution with a function 6 G D . □ 



5 Convergence of correlation functions 

Here we prove the convergence (1.8) of the correlation functions of measure /it . This implies 
the convergence of the characteristic functionals fi t in the case of Gaussian measures /i , 
H±. 

Lemma 5.1 Let S0-S2 hold. The following convergence holds as t — > oo 

Qnx,y)^QZ(x,y), Vx,yGlR 3 , Vi,j = 0,l. (5.1) 
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Proof We prove the lemma again for i = j = in the particular case, uq = almost 
surely. The general case is considered in Section 6. Let us assume for a moment that the 
function v (z) is continuous almost surely. Then the Kirchhoff formula (1.10) and the Fubini 
Theorem give 

Q^( Xi y)=Eu(x,t)u(y,t) = j^-- J dS(x') j Q^{x\y')dS{y'). (5.2) 

71 S t (x) S t (y) 

This integral is the convolution of Q]^{x,y) in both variables x, y with a distribution of 
compact support. The convolution of distributions with compact support is commutative. 
Therefore, the assumption on the a.s. continuity of v (x) can be removed by a convolution 
with a function 6 G D . Changing the variables x' = x + ut in the right hand side of (5.2), 
we get 



dS(x') J Ql\x',y')dS(y') 



(Ant) 2 

& (:r) S t (y) 

= 7^ j dS(uj) j Q l\ x + U jt,y')dS(y') + -^ J dS(u) j ' Q${x + oot,y')dS(y') 

|w|=l,W3<0 St(y) |w|=l,w 3 >0 St(y) 

= I4t,x,y) + l + (t,x,y). (5.3) 
Let us recall that v(r) = v 2 (r) . 

Definition 5.2 C U (M 3 ) is the space of functions f(y) e C(IR 3 ) s.t. \f(y)\ < Cv(\y\) with 
a constant C G 1R . 

Let us define for f(y) G ^(IR 3 ) 

nf{v) = -^- 2 J dS(co) J f(p)d 2 p, v G IR 3 . (5.4) 

|w|=l,±a) 3 >0 p-oj=v-u> 

Here d 2 p is the Lebesgue measure on the plane p ■ uj — v ■ uj . Note that the integrals with 
± are identical and converge due to (2.8). Hence, the operator 1Z : CV(IR 3 ) — > Cb(IR 3 ) is 

continuous with the obvious norm in C v : \\f\\c v — SU P ' 



yen 3 v{\y\) ' 

The convergence (5.1) follows for % — j — from (5.2), (5.3) and Lemmas 5.3 and 5.4. 
Lemma 5.3 Let S2 hold. Then for x,y G IR 3 , 

I±(t,x,y) -> 1lq± (x - y), t -> oo. (5.5) 
Lemma 5.4 Let f(y) G ^(IR 3 ) . Then 

Kf = ~£ * /. (5.6) 
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Lemma 5.4 is proved in Appendix B. 



Proof of Lemma 5.3. For a moment we assume additionally (4.4). Denote by In the 
inner integral entering (5.3): 

hi=In{x,y,uj,t)= J Ql 1 (x + ut,y , )dS(y'). (5.7) 

s t (y) 

Change the variables y' = y + out + p and denote R — \x — y \ . (4.4) implies that Q^(x + 
ut, y + out + p) = for \p\ > r + R , hence (5.7) becomes 

/n = / Ql\x + ut,y + ut + p) dS(p), (5.8) 
S t (-wt)nB 

where B denotes the ball \p\ <r + R. The sphere S t (—out) contains the point 0, hence 
in a neighborhood of the origin the sphere converges to its tangent plane u 1 - as t — > oo . 

Further, consider the case u 3 < and u 3 > separately. For w 3 < and sufficiently 
large t > t(oo) > , 

x 3 + u 3 t < -a, y 3 + u 3 t + p 3 < -a, for \p\ < r + R. 

Then SI implies that 

Q^^x + u^y + ut + p) = q l }(x-y-p) (5.9) 

Therefore, if uj 3 < , 

Jn-> | q n (x-y-p) d 2 p, t^oo (5.10) 
w-LnBo 

that coincides with the inner integral in the right hand side of (5.4), with / = q 1 } and 
v = x — y . Similarly for uj 3 > . Lemma 5.3 is proved with the additional assumption (4.4). 
At last, Lemma 4.2 and S2 give the uniform smallness of integral (5.7) over \p\ > r + R 
with large r . Therefore, (5.10) holds for any u with u 3 ^ . Hence, (5.5) follows by the 
Lebesgue Theorem on dominated convergence. □ 



6 Correlation functions in general case 

We prove Lemmas 4.1 and 5.1 in the general case. Let us assume for a moment that 
uq £ C 1 (IR 3 ) and vq £ C(1R 3 ) almost surely. Then we apply the general Kirchhoff formula 
for the solution u(x, t) to the Cauchy problem (1.1): formally, 

1 f 1 

u(x,t) = — J (y {x') + -u (x') + Vu {x') ■ n x (x'))dS(x'), (6.1) 

S t (x) 



x' — x 



where n x (x') = — . It implies similarly to (5.2), 



Q?(x,y) = -^— 2 J dS(x') J ([Q^{x\y') + V yl {V xl Q^{x\y').n x {x'))-n y {y') 



Unt) 

St(x) St(y) 
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1 

+ 1 



QSV, vO + Qo(x', y') + v„QjjV, y ') • n y ( y ') + v^QjjV, y ') • n.M + ? q°V, z/) 



+ 



V^Q°V,y') • n x (aO + V y ,Q™{x\y') ■ n y {y')]) dS(y'). (6.2) 



Proof of Lemma 4.1 in the general case We will prove the uniform bounds for Q®°(x, x) , 
V x • VyQ° t °(x,y)\ x=y and Q\ l (x,x) . Then (4.2) implies (4.1). 

Step 1 (6.2) represents Qt°(x, y) as the sum of convolutions with D^yQo l (x,y) in both 
variables x, y, with < d = k + I + \a\ + \/3\ < 2. Therefore, V x • V y Q^ (x, y) is the 
similar sum involving D"'^QQ l (x,y) with 2 < d < A. The similar representation holds for 
Q\ l (x ) y). Hence, • V y Qf ) (x ) y) and Q] 1 (x,y) can be estimated by the method of the 
proof of Lemma 4.1 in Section 4. Indeed, due to S2 with d = 2 , (2.8) and Lemma 4.2 we 
get (cf. formula (4.8)), 

V x -V y QT(x iy )\ x=y + Q\\x,x) < J v 2 {\x' -x"\)dS{x')dS{x") < d < oo. 



St(x)xS t (x) 



Step 2 Qi (x,y) requires the particular attention due to the presence in the integrand of 



the functions D^^Qq (x' , y') that are estimated by Ud(\x' — y'\) with d = 0,l. In this case 
due to (2.8) and Lemma 4.2 we have to analyze (6.2) more carefully. The corresponding 
contribution of D a f yl Q^{x\y') with d = k + I + \a\ + = 0, 1 is 



1 



Q™( X ',y') + ... + -Q™( X ',y') 



dS(y'). 



dS(x') J 

St(x) S t (y) 

Lemma 6.1 The integral I^°(x,y) converges to zero as t — > oo . 
Proof The assumption S2 implies 

l 7 °V,y)l<7^L / dS(x') J - t [AM\x'-y'\) + - t MW-y'\)]dS(y'). 



Unt) 

St(x) 

Therefore, Lemma 4.2 implies 



(6.3) 



St(y) 



C c 1 r r 2t I 

I™(x,y)\ < — J dS(x')- t [J o (ArMr) + - t rMr))d 



St(x) 



Now (2.8) implies the convergence to zero by the Lebesgue theorem. 
Lemma 4.1 is proved in the general case. 



(6.4) 

□ 
□ 



Proof of Lemma 5.1 in the general case We will consider i — j — . The other cases 
can be considered similarly. 

\P\ < 1 entering (6.2), 



Step 1 The integrals of ,Qq 1 (x / ,y') with d = k + I 



\a\ 



converge to zero by Lemma 6.1. For the integrals of D"r yl Q^ l (x' , y') with 2 < d < 4 , the 
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convergence follows by the method of the proof of Lemma 5.3. Let us define for the functions 
/ G C*(IR 3 ) := {/ G Ll c (M 3 ) : |V/(y)| G C^M 3 )} (cf. Definition 5.2), the operator 

Vf{v):= T^ J dS{i0) J W(P)-^^, ^GIR 3 . (6.5) 

M=l,w 3 >0 v-u>=p-u> 

With the obvious norm in C x v : ||/|| C i = sup ^.{^ , the operator V : C3(IR 3 ) -> C 6 (IR 3 ) 

" s/eR 3 KM) 

is continuous. For instance, the operator "P can be applied to q± with 1 < k + 1 < 2 since 
G C^(IR 3 ) by S2. Similarly, the operator 1Z (see formula (5.4)) can be applied to D a q± 
with 2 < k + I + \a\ < A since G C^IR 3 ) . Now, (6.2) and the method of proof of 

Lemma 5.3 imply the convergence (5.1) with % — j — to the limiting function 

g f = n[q? + q n - A(g°° + q 00 )] +v[q° + 1 - q m - q% + q]_ ] . (6.6) 

Step 2 It remains to prove that = q^ . First, let us prove that 

ttAtf = -\<% ■ (6-7) 



1 



In fact, A(7£Ag°°) = — -Ag°° due to (5.6), hence f(x) = IZAq ® — q Q + is a smooth harmonic 

function in IR 3 . On the other hand, Aqf G <^,(IR 3 ) by S2. Hence, g(x) = IZAqf G 
Cb(JR 3 ) , and moreover, 

g(x) — > 0, |x| — > oo. (6.8) 

Indeed, 

oo 

I J Aqf{p)d 2 p\< J v{\p\)d 2 p = 2tt J ru(r)dr (6.9) 

p-U) = X-LO p-U>=X-UJ X-U) 

similar to (4.7) with t — oo . This integral is bounded uniformly in \u\ — 1 and converges 
to zero if \x\ — > oo and x = \x\9 with 9 ■ uo ^ . Therefore, (6.8) follows from (5.4) by the 
Lebesgue theorem. Further, \f(x)\ < \g(x)\ + ^o(M) again by S2. At last, vo(r n ) — > for 
some sequence r n — > oo due to (2.8). Finally, the maximum principle and (6.8) imply for 
any fixed x G IR 3 , 

\f(x)\ < max \g(y)\ + z/ (r n ) -> 0, n -> oo. 

|j/|=»Vi 

Therefore, /(it) = and (6.7) is proved. Further, let us consider the terms with V in (6.6). 
Obviously, is a convolution. We prove the next lemma in Appendix A. Let us recall 

thatP(oO = -^ SgnA;3 - 



\k\ 

Lemma 6.2 For f G D we have 

Vf = \p*f. (6.10) 

Let us assume for a moment that all the correlation functions q± (•) are smooth and have 
a rapid decay. Then (6.6) coincides with (2.13) by (6.7) and Lemmas 5.4, 6.2. In the general 
case we consider the formula (6.10) as the definition of the convolutions with P , entering 
(2.13)-(2.15). Lemma 5.1 is proved in the general case. □ 
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7 Bernstein's argument for the wave equation 



In this and the subsequent section we develop a version of the Bernstein 'room-corridor' 
method. We use the standard integral representation for the solutions, divide the domain 
of integration into 'rooms' and 'corridors' and evaluate their contribution. As the result, 
(U(t)Y , ty) is represented as the sum of weakly dependent random variables. We evaluate 
the variances of these random variables that will be important in next section. 

For the wave equation the similar method has been used in [10, Section 6] for an odd 
n > 3. Our mixing condition S3 is different from [10] (see Remark 2.7). Respectively, the 
method of [10] requires a suitable modification. 

Denote by £t(x) = £(x,t) = — <KM 2 ~ t 2 ) the fundamental solution to the wave equa- 

2% 

tion. The support of St is the sphere St = {x e IR 3 : \x\ = t} . Therefore, the dynamical 
group U(t) of the problem (1.2) is the convolution operator 

U(t)Y = Qt * Y , t > 0, (7.1) 

where 

( E t £ t \ (7 2) 



Next we introduce a 'room-corridor' partition of the space IR 3 . Given t > , choose 
d = d t > 1 and p = p t > and an integer N = N t > . Asymptotic relations between t , 
d t and p t are specified below. Define 

ai = —t, bi — ai + d; a 2 = b x + p, b 2 = a 2 + d; . . . , b N = a N + d = t. (7.3) 

We divide the sphere S t by the planes orthogonal to the axis Ox 3 into the slabs which we 
call the "rooms" R\, (k = 1, N) , separated by the "corridors" C\ (k = 1, N — 1) , 

R{ = {xeS t : x 3 E [a k , b k }}, Ci = {xeS t : x 3 e [b k , a k+1 }}. (7.4) 

Here x = (xi,X2,x 3 ) , d is the width of a room, and p of a corridor. Then 

St=(uBt)u(\JCf). (7.5) 
For any region E C S t we define the distribution £ t s with the support in E 

(St^6):=^-j6(z)dS(z), 6eD. 
Note that Vt > £ t (x) = --S\\z\ 2 - t 2 ) = -E t (x) - V(-S t (x)). For any region EcS t 

7T t t 

we define the distribution £\y. '■ 

E^{x)^\s t ^ X )-V{- t S t ^x)). (7.6) 



Then for E = S t we have £ t ,T, = £t ■ Let us denote 



^ := ^ • (7.7) 
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We define the random variable 

/ t (E) = (&,E*>o,*>, (7.8) 
where $ 6D is a fixed function from (2.17). For instance, define 

r k ^I t (R k ), c k = I t (C k ). (7.9) 

(7.5) implies that 

N t N t -1 

(U(t)Y ,9) = (gt + Yo,*) = £rf + E^- (7.10) 

k=l k=l 

Lemma 7.1 Let SO, S3 ZioZrf. The following bounds hold for t > 1 and V/c 

£|r*| 2 < C(*H/t, (7.11) 
£|c^| 2 < C(tt) a/1 (7.12) 

Proof. We prove the following estimate: for any region ScS| 

E\I t {^)\ 2 < C(^)|E|/t 2 . (7.13) 

Then (7.11) and (7.12) would follow from this estimate with E = R k and E = R k , respec- 
tively, as \R k \ — 2ntd t and \C k \ = 2ntp t . 

Now we prove (7.13). From (7.7) and (7.8) it follows that for * = ^ x ) e V 

/(E) = (4 E * uq, *°) + (£ t , s * «o, *°> - (4s * Vm , W 1 ) + (4s * v , V 1 ). (7.14) 
Substituting (7.6) in the first and the last terms in the RHS of (7.14), we get 

/(E) = (^ s *v^* )-((f^s)*v«o,* )-(^,E*v«o,v* 1 ) 

+(£*,£ * vo, *°) + -(4e * ^o, + ((f 4e) * ^o, W 1 ). 

Hence, 

M 

/ t (E) = ]T //, where // = c,-(f)(^ E * 0,-). (7.15) 
j'=i 

Here M < 6 , Cj(t) is a bounded function for t > 5 > , Wj is one of t~ x uo , V«o or t> , 
4e is one °f 4s C 2 ) or -£t,v{z) , % is one °f D a ^/ 1 with |a| < 1 . Therefore, for t > 1 

M M _ 

3=1 3=1 

M l r r,\ 



< CiEsK/ ^o°(^-^-(z/-p))(7.i6) 

i=l S E l«l = l/3|=l 

+g 1 (x-z -( y - P ))) ds(z)ds(p), e,{x)e 3 {y))\. 
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We have 

supp * C B ro = {x G H 3 : \x\ < r } (7.17) 

with an r > . Since x, y G supp 0j C supp ^ C £> ro , \x — z — y + p\ > (\z — p\ — 2r )+ , 
where s + = max(s, 0) , s G IR. Since z/ d (r) are non- increasing functions, (7.16) and S3 
imply 

£|/ t (E)| 2 < C(*)i / dS(z) J (^u ((\z - P\ ~ 2r ) + ) + ~ p\ ~ 2r )+)) dS(p).(7.18) 

s s 

Then Lemma 4.2 and (2.8) imply as in Lemma 6.1, 

. . 2i 

£|/ t (£)| 2 < cU | (ri*((r - 2r ) + ) + £n,((r - 2r )+)) rfr < . □ 



8 Convergence of characteristic functionals 

In this section we complete the proof of Proposition 2.10. If Qoo{^, ^) = Proposition 
2.10 is obvious, due to (5.1). Thus, we may assume that 

Q oo ^,^)^0. (8.1) 

Choose < 5 < 1 , and 

^-(ln^+l)) 1 / 10 , p^t 1 -*, t^oo. (8.2) 
Lemma 8.1 The following limit holds true: 

N t (v 2 (p t ) + (^) 1/2 )+N?(v (p t ) + ^)^0, t^oo. (8.3) 

Proof. Since u d (r) are non- increasing functions, (2.8) implies 

Mr)Hr + l)=[^ds< [^±ds<C <oo. 



Then (8.2) implies (8.3). □ 
By the triangle inequality, 

|&(tt)-/U*)l < |^exp{i(C/(0y o ,*>}-^exp{^ t r t fe }| + 

+ l ex P{-^E t ^(^) 2 } " exp{-iQoo(^,^)}| + 

+ |£exp{^rf} - ex P {-i£^(f ) 2 )l 
= I 1 + I 2 + I 3 , (8.4) 
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N t 

where the sum J2t stands for J2 ■ We are going to show that all the summands I\ , h , ^3 

k=i 

tend to zero as t — > 00 . 
Step (i) Eqn (7.10) implies 

h = |£exp{i£ t r t fe }(exp{^c t fe } - 1)| < ^E\^\ < £ t (£|cf| 2 ) 1/2 . (8-5) 

/From (8.5), (7.12) and (8.2) we obtain that 

h < CN t (p t /tf 2 -> 0, t^oo. (8.6) 

S'tep ^ By the triangle inequality, 

h < ^\Y,E(rt) 2 -Quy,n<l\Qt(y,y)-Qoo^,m 

= / 2 l+/ 22 + /23, (8.7) 

where Q t is the quadratic form with the integral kernel (Q l t( x i y)^j . Eqn (5.1) implies 
hi — > 0, t — > 00 . As to J 22 , we first obtain that 

^ = ^ KE/f) 2 - E^) 2 I < E I^MI- (s-s) 

The next lemma is a corollary of ([11, Lemma 17.2.3]). 

Lemma 8.2 Let £ be a random value measurable with respect to the a -algebra a dl {A) , i] 
be a random value measurable with respect to the a -algebra o- d , 2 {B) , and dist(A, B) > h > . 
1) Let {E\i\ 2 fl 2 < a , {E^ 2 ) 1 ' 2 <b . Then 

\E£ V - EiE V \ < C ab (ff A (h). 

ii) Let |£| <a, \q\ < b almost surely. Then 

\E£ri - EZE V \ < C ab <f> dl , d2 (h)- 

We apply Lemma 8.2 to deduce that J 22 — > as t — > 00 . Note that r\ = (G tjR k *Y , 
is measurable with respect to the a -algebra a dk (A k ) , where 

A k = {x - y : y E R k , x E supp * c B ro }. 

The distance between the different rooms R k is greater or equal to p t according to (7.3) 
and (7.4). Then p(A k ,A l ) > p{R k ,R\) - 2r > p t - 2r . Hence (8.8) and SO, S3 imply, 
together with Lemma 8.2 i), 

hi < CN 2 u ((p t - 2r )+) -> 0, t^oo, (8.9) 
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because of (7.11) and Lemma 8.1. Finally, it remains to check that ^23 - ► , t — > oo . By 
Cauchy-Schwartz inequality, 

'» < l^(E^) 2 -^(E^ + E^) 2 l 

< N t ^E\c^ + 2(E(^) 2 ) 1/ \N t ^E\c^) 1/ \ (8.10) 

(7.11), (8.8) and (8.9) imply E(J2 t r k ) 2 < C 1 + C 2 N?v ((p t -2r ) + ) < C 3 < oo . Then (7.12), 
(8.10) and Lemma 8.1 imply 

/ 2 3<CW t V* + ^V t W0 1/2 ^0, t^oo. (8.11) 

So, 

I21 , -^22 , -^23 tend to zero, as t — > oo . Then (8.7) implies 

/ 2 < \ lE*S(r?) 2 - Qoo(^, *)| - 0, t - oo. (8.12) 
ifftep (raj It remains to verify 

7 3 = |£;exp{^r t fe }-exp{-^£;(r t fc ) 2 }H0, f - oo. (8.13) 
Using Lemma 8.2, ii) we obtain: 

Aft 

\EeM^ t r t}-UEeM^}\ 
k=l 

N t N t 

< \Eexp{irl} exp{i ]T r k } - E exp{ir]}E exp{i ^ r k t }\ 

k=2 k=2 

Nt N t 

+ \Eexp{irl}Eexp{iY^r k t }-\{Eexp{iT k t }\ 

k=2 k=l 

Nt N t 

< vl{{ Pt - 2r ) + ) + \Eexp{t^r k } - J[ Eexp{ir k }\. 

k=2 k=2 

We then apply Lemma 8.2, ii) recursively and get, according to Lemma 8.1, 

N t 

\Eexp{{£r k t ] - J] Eexp{ir k }\ < N t u 2 ((p t - 2r )+) - 0, t - oo. (8.14) 

k=l 

It remains to verify the convergence 

in^exp{^}-exp{--^E(rf) 2 }|^0, t ^ oo. (8.15) 
k=i 1 

According to the standard statement of the Central Limit Theorem (see, e.g. [17, Thm 4.7]) 
it suffices to verify the Lindeberg condition: We > 

-E t ^KT->o, ( 8 - 16 ) 
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Here a t = J2tE( r l) 2 > an d -£5/ = E(Xgf) , where X$ is the indicator of the event |/| > 5 2 . 
Note that (8.12) and (8.1) imply 

<7 t -> Qoo(^,^)^0, t^OO. 

Hence it remains to verify that Ve > 

£ t £*KT->0> t^oo. (8.17) 

We check (8.17) in Sections 9, 10. Finally, (8.4) and (8.6), (8.12)-(8.15) imply Proposition 
2.10. □ 



9 The Lindeberg condition 

The proof of (8.17) can be reduced to the case when for some b > we have, almost surely 
that 

|lo(aOI<&, xeJR 3 . (9.1) 

The general case can be covered by the standard cutoff argument in the following way. 
We decompose Y in two summands: the first one, satisfying the estimate (9.1), and the 
remainder. For large b , the dispersion of the remainder is small due to S2, S3 and Lemma 
8.2, i), then the dispersion (7.11) of the corresponding variables r\ is small uniformly in t . 
The last fact follows from the proof of (7.11). 
Further, we estimate 

Y E £ \r k t \ 2 = Y litfl JL-EJr*] 2 < 4nt 2 max -\-E £ \r k t \ 2 . (9.2) 
Therefore, it remains to prove 

max -^-EJrfl 2 = o(r 2 ), t -> oo. (9.3) 
k=i,..,N t \R*\ £l tl v ; ' K J 



The Chebyshev inequality implies 



E £ \r k t \ 2 < —E\r k \ 4 . (9.4) 



Using (7.15), we get 

E\r k t \ A = E\I] + ... + J t M | 4 < C{M)E(\ll\ A + ... + |/f| 4 ). (9.5) 

Therefore, (9.3) follows from the estimate 

1 



max Jj^ E \ < £t,R$ *w k ,0 k >\=o(t ), t -> oo. (9.6) 
We prove the following proposition in the next section. 
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Proposition 9.1 Let (9.1) holds, and w = t 1 u ,Vu or v . Then for any S C St the 

bound holds 

E\<£ t ^*w,e>\ 4 <C(6)(jJ |S| 2 . (9.7) 

— ■ Z 

Here E 3 tTi is one of £t,z(z) or —£t&(z) j Oj is one of D a ^ 1 with \a\ < 1 . 

t 

This proposition implies (9.6): 



^E\ < £ t ^*w k ,6 k > | 4 < j^cw (j) \R$\ 2 < c(&,*)!31 



since \R\\ < Ant 2 /N t , where N t — > oo . (8.17) is proved. □ 

10 The fourth order moment functions 

We deduce Proposition 9.1 from the bounds for the fourth order moment functions. 

Denote by mf(z) := Ew(zi) ■ ... ■ w(zi) , z = (zi, z{) , where w(zk) = Vo(zk) for every 
k = I,..., I, or w(zk) = Vuo(zk) for every k = I,..., I, or w(zk) = t^u^Zk) for every 
k — 1, I . We have supp# C B ro for an r > . Then left hand side of (9.7) is estimated 
as follows, 

E\ <Stp*w,0> | 4 < j j \m ( o\x-z)\dS(z)dx, (10.1) 

where dS(z) := dS(z\) . . . dS{zi) . Therefore, we have to prove that 

I(x) = J \m { o\x-z)\dS(z) < Cb 4 m 2 , xeB ro . (10.2) 

Step 1 Let us prove an estimate for the moment functions m^\yi, y 2 , 2/3, 2/4) by the method 
[10, Section 6.2]. We use the mixing condition for different configurations of the points 
2/1,2/2,2/3, 2/4 in the space IR 3 . 

Lemma 10.1 The bound holds 

|4 4) (2/i,2/2, 2/3, 2/4)1 < 46 4 (v 2 2 {\\ yi -2/2I) +^l2/i -y 2 |)r 4 ) (10.3) 

+ 1664 E (^2-iflj/l - 2/3 1) • - J/4 1) + ^2 2 - l (l2/l - 2/4|) • ^2 2 -i(l2/2 ~ 3/3 1 

0,2 

Proof. Let us divide the space H 3 in three regions I±, h by two hyperplanes that are 
orthogonal to the segment [2/1,2/2] and divide it in three equal segments, y\ e I\ , 2/2 £ ^3 • 
At least one of the regions ii, 7 2 , ^3 does not contain y 3 , y 4 . If the points y 3 , 2/4 ^ ii , then 
SO and S3 imply (10.3), since 

|m 4) (2/i, 2/2,2/3,2/4) I = |mo 4) (2/i,2/2,2/3,2/4) -mf'd/irf (2/2, 2/3, 2/4) I 
< 46 4 (i^(i|y 1 -i te |) + i^(i|y 1 -y i! |)r 4 ). 
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The same proof is valid for the case 1/3,1/4 G" h . Now let us assume that 1/3,1/4 G" h , for 
instance, 2/3 G h , 1/4 & I3 . Then SO, S3 imply (10.3), since by Lemma 8.2, ii) 

I racftyi, 2/2,2/3,2/4)1 < l"^o 4) 2/2, 2/3, 2/4) -"io 2) (2/i,2/3)^o 2) (2/2, 2/4)1 + |^ 2> (2/1 , 2/s)^ 2) (2/2, 2/4) I 
< 4& 4 (vl{\\Vi - 2/2I) + ^(|l2/i - 2/ 2 |)t- 4 ) + 166 4 £ ^ ^-^(bi - 2/ 3 |) • ^ 2 %(|2/2 - y*\). 
The proof for the case 1/3 E I3, 2/4 G ii is the same. □ 

Remark For a translation-invariant measure /x the estimate similar to (10.3) is obtained 
in [1, inequality (20.42)]. 

Step 2 (10.3) holds with any permutations of 2/1,2/2,2/3,2/4 i n the RHS. Hence 

|4%)| < 4&4 (>4(\\v. -y P \) + ^o(|l2/, - y P \)t- 4 ) 



i,j=0,2 
rl 1 i\/f2 /- 



+ 16&4 S {Ti u 2-i(\ys - yk\) ■ T-A-j{\y P - yi\) + j^i-idy* - yi\) ■ ^Lj(\y P - Vk\ 



= Ml p (y) + M^y) (10.4) 
for any permutation {s,p, k,l} of {1,2,3,4}. Let us define 

S s ^p i = |^ G S J J z s Zp J = max | z% Zj | } . 



Then (S) 4 = (J S S)P , where the union is taken over all the pairs (s,p) of the indexes 

(*.p) 

1,2,3,4. Therefore, (10.4) implies 
J(x) = /|mj 4) (x-z)|^)<^{ I Mj^x - z)dS(z) + I M s 2 p (x-^(z)}. (10.5) 

S 4 ( S 'P) s s 



Here the sum is taken over all the pairs (s,p) . Every of the six terms corresponding to 
different pairs (s,p) in the RHS of (10.5) coincide. We have to estimate I(x) only for 
x G -B 4 (see (10.2)). Then \z s — z p — x s + x p \ > (\z s — z p \ — 2r )+ for any z S) z p G IR 3 . Since 
u d is a non-increasing function, (10.4), (10.5) imply 

< Cb' J - z 2 \ - 2r )+) + ui(^(\ Zl - z 2 \ - 2r ) + )r 4 ) dS(z) 

Sl,2 

+Cb 4 g / " *l - 2r ) + ) • ^ 2 %((k - *| - 2r )+)) dS(*) 

*>j— 0,2 Si2 

= Ji + J 2 . (10.6) 
S'tep 3 Let us estimate Ji and J 2 separately. 
Lemma 10.2 h < C6 4 |S| 2 . 
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Proof The integrand in I\ does not depend on z% and z^ . Therefore, the result of the 
integration in the Z3, Z4 we estimate by the factor 7r|E||zi — 2:2 1 2 , since \z^ — 24 1 < \zi — z 2 \ 
by the definition of Ei^ • Lemma 4.2 implies 

h < Ctfm /" - ^ 2 | - 2r ) + ) + - - 2r ) + )) |^ - ^ 2 | 2 rfS'C^OrfS'C^) 



2 



< Ci6 4 |E| 2 J (v 2 2 {\{r - 2r )+) + rt 2 (^(r - 2r )+)) r 3 dr. (10.7) 


(2.8) implies 

r r 

r 2 z/ 2 (r) = ^ 2 ( r )2 y sc?s < 2 J sv 2 (s)ds < C < 00. 


Therefore, using (2.8) again, 

2t 

J r 3 v 2 (-(r -2r )+)dr < C < 00. (10.8) 


Finally, the integral 

2 / , 1 1 , **bd(r-2r ) + ) n,(I( r -2r ) + )r 4 
y u 2 (-(r-2r ) + )-r 3 dr = J —3— 3 _ rfr 



is bounded: it converges to zero as t — > 00 by Lebesgue theorem as in Lemma 6.1. Hence, 
Lemma 10.2 follows from (10.7) and (10.8). □ 

Lemma 10.3 I 2 < C6 4 |E| 2 . 

Proof Since v 2 d {j) < Cv d {r) , and E 1)2 C E 4 , we have by Lemma 4.2 
/ 2 = C6 4 E y ^2-i(ki - z 3 \ - 2r ) + ) • \ 3 v 2 _ 3 {\z 2 - z 4 \ - 2r ) + ) dS(z) 
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< C6 4 ^ y ^ i/ 2 -i ( I ^1 — ^3 1 — 2r ) +) diS' d5 (^ 3 ) y ^-v 2 _ j (\z 2 -z 4 \-2r ) + ) dS^dS^) 

i,j=0,2 S 2 S 2 
2t 2t 

< Cb* ]T l S l y ^2-*((r - 2r ) + ) rfr • |E| J ^2-j((r - 2r )+) dr < C5 4 |E| 2 . (10.9) 

i,j=0,2 o 

In the last inequality we use (2.8) and the Lebesgue theorem as in Lemma 6.1. Lemma 10.3 
is proved. □ 

Now Lemmas 10.2, 10.3 and (10.6) imply (10.2). □ 
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11 Convergence to equilibrium for variable coefficients 

We extend all results of previous sections to the case of the wave equations with variable 
coefficients. We consider the wave equations in IR 3 with the initial conditions 

3 

u(x, t) — E dj(a jk (x)d k u(x, t)) - a (x) u(x, t), x G IR 3 , t G M, 

j,k=i (11.1) 

u\ t=0 = u (x), u\ t =o = v (x), 
d 

where <9, = - — . We assume the following properties E1-E3 of Eqn (11.1). 

OXj 

El ajk(x) = 5jt + bjk(x) , where bjk(x) G D ; also a (x) G D . 
E2 a (x) > 0, and the hyperbolicity condition holds: 3a > s.t. 

1 3 

H(x, k) = -J2 Oij(x)kikj > a\k\ 2 , x, k G IR 3 . (11.2) 

E3 Non-trapping condition holds, [22]: for (x(0), fc(0)) G IR 3 x IR 3 with fc(0) ^ 

|x(t)| — > oo, t — > oo, (11-3) 

where (x(t), /c(t)) is a solution to the following Hamiltonian system 

x(t) = V k H(x(t), k(t)), k(t) = -V x H(x(t), k(t)). 

Example. E1-E3 hold in the case of constant coefficients, a,jk{x) = 5ij . For instance, E3 
hold because k(t) = =^ x(t) = fc(0)i + x(0) . 

We denote as above, F(t) = t), t)) , Fo = ( M o,^o) • Then (11.1) becomes 

Y{t) = ^{Y{t)), t g ir, y(o) = y - (ii.4) 

Proposition 2.2 holds for the solutions to the Cauchy problem (11.4) as well as for (1.2). 
Let Yq in (11.4) be a measurable random function with values in (Tl, B(7i)) , and let Hq be 
its distribution, as above. Denote by fi t the distribution of the solution Y(t) to the problem 
(11.4). Let us state the extension of main Theorem 2.8. We introduce the appropriate Hilbert 
spaces of initial data of the infinite energy. Let 5 be an arbitrary positive number. 

Definition 11.1 H$ is the Hilbert space of functions Y — (u,v) G Tl with the finite norm 
|||Y|||2 = J e' 2S ^(\u(x)\ 2 + \Vu(x)\ 2 + \v(x)\ 2 ) dx < oo. 

Theorem 11.2 Let E1-E3, S0-S3 hold. Then 

i) the convergence (2.4) holds for any e > . 

ii) The limit measure fi^ is a Gaussian measure on 7i . 
Hi) The limit characteristic functional has the form 

/i 00 (^) = exp{-iQ 00 (^vI/,^xI/)}, y eV , 

where W : T> — > 7i' s is a linear continuous operator for sufficiently small 5 > . 

Theorem 11.2 follows immediately from Theorem 2.8, using the method [10]. The method 
is based on the scattering theory for the solutions of infinite energy. 
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12 Appendix A. Radon transform 

Proof of Lemma 5.4 Since J f(p)d 2 p is an even function with respect to u , it suffices 



p-LU = Z-LU 

to prove the next lemma. 

Lemma 12.1 Let (2.8) hold, and f G ^(H 3 ) . Then 

j±- 2 j dS{u) J f( p )d"p = - 1 -S*f(z), V^GlR 3 . (12.1) 

|w|=l p-LU=Z-LU 

Proof. Both sides of (12.1) define the continuous operators C U (M 3 ) i— > Cf,(IR 3 ) . Therefore, 
it suffices to consider / G D . Applying the Fourier transform, we obtain with p = \k\ , 

(£ *f)( z ) = —J £(k)f(k)e^ k d 3 k = / d5( W ) / pV^(pc,)/(p^P- 

|w|=l 

(12.2) 

1 

We substitute S(pu) = — - in the right hand side of (12.2) and get 

+oo 

(S*f)(z) = -—- J dS(uj) J e-^f{pu)dp. (12.3) 

^ n) \u\=l 

Note that 

+oo 

f(pu)= j e^ h f(h,u)dh, where f\h,u)= J f(y)d 2 y. (12.4) 

— oo y-u)=h 

Then from (12.3), (12.4) we have 

+oo +oo 

(£*/)(*) = -_ / d5(w)- / e -*- w dy / e^ h f(h,uj)dh 



|«|=i 



= -g^ / ^H^mW^M) = -g^ / /"(*-a;, W )dS(a;). 

|w|=l |«|=1 

Lemma 12.1 is proved. □ 
Proof of Lemma 6.2. Since F[P](/c) = — — sgn/c 3 , we have 

+00 

(P*f)(z)^^JP(k)f\k)e- iz - k d 3 k^--^ 3 f dS(u) J pe-^-sgnM/Mdp 

|w|=l 
+oo . +oo 

= -T^p / dSM / pe-^/Mdp+^s / / Pe- ipz0J f(f^)dp. (12.5) 

|w|=l,w 3 >0 |w|=l,w 3 <0 
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In the last integral we change the variables u — > — uj , p — > — p , then apply (12.4) and get 

+00 

(P * /)(*) = -tAs / / Pe^fMdp 

|w|=l,W3>0 -OO 



+OO +OO 



J dS{uj) J e~ ipz -" pdp J e iph f(h,u)dh. (12.6) 



(2tt) 3 

|w|=l,W3>0 -OO -OO 

Note that 

+00 +00 

P J e iph f(h, u)dh = i J e iph {Vf)\h, u)-udh, p G IR. (12.7) 
—00 —00 
Indeed, applying (12.4) in the both sides of 

F[Vf](pu).u = -ipF\f](pu), 

we obtain (12.7). Finally, from (12.7) and (12.6) we get 

+00 +00 

(P*f)(z) = 7^3 / dS{u) J e-^dp J e ph (yff(h,u)-udh 



|w|=l,W3>0 -OO 



|a)|=l,W 3 >0 

= 4^2 / (V/)*(z-a;,a;)-a;dSM 



|w|=l,a) 3 >0 

= 4P/(z). (12.8) 
Lemma 6.2 is proved. □ 

13 Appendix B. Gaussian measures in Sobolev's spaces 

We verify (3.3). Definition (3.2) implies for u G H s a , 

|| u ||2 q = J{1 + \x\) 2a ( J e-^-^l + |£;|) s (l + \n\) 8 u(k)ti(ri) dkd V )dx. (13.1) 

Let p(du) be a Gaussian translation invariant measure in H s a with a correlation function 
Q(x,y) = q(x — y) . Let us introduce the following correlation function 



C(k,rj) = J u(k)u(rj)n(du) (13.2) 

in the sense of distributions. Since u(x) is real- valued, we get 

C(k, 77) = F x ^ k F y ^ v Q(x, y) = CJ(k - V )q(k). (13.3) 
Then, integrating (13.1) with respect to the measure p(du) , we get the formula 

/ \\u\\l a p{du) = C n J(l + \x\) 2a dxJ(l + \k\) 2s q(k)dk. (13.4) 

Substituting q(k) = 1 and q(k) = \k\~ 2 , we get (3.3). □ 
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